Real-Time Ginzburg-Landau Theory for Bosons in Optical Lattices 
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Within the Schwinger-Keldysh formalism we derive a Ginzburg-Landau theory for the Bose-Hubbard 
model which describes the real-time dynamics of the complex order parameter field. Analyzing the ex- 
citations in the vicinity of the quantum phase transitions it turns out that particle/hole dispersions in the 
Mott phase map continuously onto corresponding amplitude/phase excitations in the superfluid phase. 
Furthermore, in the superfluid phase we find a sound mode, which is in accordance with recent Bragg 
spectroscopy measurements in the Bogoliubov regime, as well as an additional gapped mode, which seems 
to have been detected via lattice modulation. 



Within the last decade ultracold atoms in optical lattices 
[01 0] have become a standard tool for studying quantum- 
statistical many-body effects. Due to their high tunability, 
these systems represent an almost perfect test ground for a 
large variety of solid-state models. In particular, the experi- 
mental observation of the seminal quantum phase transition 
from the Mott insulating (MI) to the superfluid (SF) phase, 
exhibited by a single -band Bose-Hubbard (BH) system of 
spinless or spin-polarized bosons, has recently attracted a 
lot of attention [3]. Although the occurrence of this quan- 
tum phase transition is evident from the momentum dis- 
tributions of destructive time-of-flight measurements, its 
precise location cannot be determined from them. Re- 
cently, however, more detailed information about the col- 
lective excitations of this system could also be achieved 
in a non-destructive way by exciting the system via lat- 
tice modulation [4] or by Bragg spectroscopy JB)]. Deep in 
the SF phase, the observed gapless excitation spectrum can 
be well understood within a Bogoliubov theory [6]. Ap- 
proaching the phase boundary, both a slave -boson method 
13] and a random-phase approximation [8] have even pre- 
dicted an additional SF gapped mode. When the MI phase 
is reached, both SF modes turn continuously into particle 
and hole excitations which are also found by mean-field 
theory [6]. By being both gapped, they characterize the in- 
sulating phase in a unique way. Until now, however, there 
exists no theory which describes the collective excitations 
in these physically different regimes in a unified way ]9|]. 
To this end we derive in this Letter a real-time Ginzburg- 
Landau theory within the Schwinger-Keldysh formalism. 
Following Refs. our approach is technically based 

on resumming a perturbative hopping expansion for the 
effective action. This is physically justified as a dimen- 
sional rescaling of the hopping parameter J, which reads 
J -> J/D for dimensions D > 1 fiHl in the presence of 
a condensate, turns out to suppress all nth order hopping 
loops by a factor l/.D ra_1 . Resumming the 1 -particle irre- 
ducible contributions up to the nth hopping order yields an 
effective (l/.D)-expansionuptothe (n-l)th order. By re- 



stricting ourselves in this Letter to the lowest order n = 1, 
we will get an approximation which is exact for infinite di- 
mensions or infinite-range hopping lfl3ll and turns out to 
describe well the quantum phase transition for D = 2,3. 
Although our approach is originally designed for describ- 
ing the vicinity of the phase transition, it turns out to repro- 
duce also the Bogoliubov spectrum deep in the SF phase. 
This is an evidence that our Ginzburg-Landau theory can be 
applied for quantitative predictions within a broad range of 
system parameters. 

For our perturbative calculation we use the Dirac picture, 
where the dynamics of the operators is given by an unper- 
turbed Hamiltonian H . To this end we introduce the split- 
ting of the BH-Hamiltonian H B h = H + H kin , where 



2 a i a z 



is the local in- 



H = E iL _ 

teraction and H kin = — ■ J^-aja^ the hopping term. 

Here a, (at) denote the bosonic annihilation (creation) op- 
erators at lattice site i, n the chemical potential, U the 
on-site interaction parameter, and the hopping ma- 
trix element being equal to J > for nearest neighbors 
only. In order to break the underlying U(l) symmetry 
of the BH-Hamiltonian, we add an artificial source term 

Hs(t) = J2i 3i(t)^i + c - c - to me Hamiltonian with ex- 
ternal currents ji(t). Considering the source term as part 
of the perturbation Hi(i) = H kin + H s (£), we arrive at 
H(i) = H + Hi(t). 

Within the Schwinger-Keldysh formalism the Dirac picture 
involves a time evolution along a closed real-time contour 
which starts and ends at some initial time t and must ex- 
tend to any large time which can be chosen to be +oo. To 
define the position of the operators on the contour, we must 
provide them with an additional path index P = ± Ill4ll . 
The resulting time-evolution operator along the closed con- 
tour then reads 



S' S = T c exp 



5>" 



dt' H 1± (t') 



(1) 



2 



where the contour-ordering operator T c turns the opera- 
tors of the forward path time-ordered, followed by the anti- 
time-ordered operators of the backward branch. If we dis- 
tinguish accordingly the currents j on the forward path of 
the contour from the backward ones, the resulting generat- 
ing functional Z[j,j*] = (t c {S^S}^ defines the corre- 
sponding Green's functions. Denoting the set of variables 
{i, t, P} by Greek indices, the contour-ordered Green's 
functions read: 



G 



■n+m— 1 



(2) 



The angle brackets (-) denote a thermal average with re- 
spect to the unperturbed Hamiltonian H . In order to ac- 
count for perturbative contributions to the thermal average, 
one should actually choose a time-contour which also con- 
sists of an imaginary part from t to t — ih{3. However, it is 
widely believed in the literature that, if we push t — > — oo, 
this imaginary part can be neglected for initially uncorre- 



cted systems II15L llfll . We will follow this tradition and 
discuss its limitations at the end of this Letter. 
The coefficients of an expansion of the generating func- 
tional Z with respect to both the currents and the hop- 
ping parameter are given in terms of unperturbed Green's 
functions. A simple diagrammatic construction rule for 
the expansion is found, when we consider the functional 
J-[j,j*] = — ilnZ\j, j*]. According to the linked-cluster 
the expansion coefficients for this functional are the cumu- 
lants or connected Green's functions [fTvL [l8h . which are 
related to the Green's functions (0 by decomposition for- 
mulas. Defining j a = P a j iaPa (t a )/h, where the sign P Q 
takes into account the direction of the time evolution, and 



in J: 



Ji a i„ j h, we get in fourth order in j and first order 



J 7 — ja [Cq/3 + P 'k J 'apC ' aK C ' K fi\ jp + jaj/3 
— CajSjS + Pfc — {C aKll sCf3 K + h.c. 



(3) 
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Doubly occuring indices must be summed (integrated). 
The above notation obscures the locality of the unperturbed 
cumulants C in its spatial variables. The general structure 
of Eq. (01 does not change in Fourier space with the vari- 
ables {k, Lo, P}, if we define J a p = b^^^pJ jh. Note that 
the frequency is conserved at each cumulant, when trans- 
formed in Fourier space. 

The relation H5J- /5j^ = (& a ), which defines the order 
parameter field motivates to perform a Legendre trans- 
formation: 



(j Q *; + c. c .) 



(4) 



In order to calculate the effective action T as a power se- 
ries in J, ^ and we must invert ^[j, j*] iteratively 
in both the order parameter fields and the hopping. For 



describing the symmetry-broken behavior, the Ginzburg- 
Landau theory demands a fourth-order term in *S>. Already 
a first-order expansion in the hopping yields mean-field re- 
sults, since all diagrams without lattice loops turn out to 
be resummed due to the Legendre transformation. This 
can be seen by noticing that the second-order term in Eq. 
d3| is local in Fourier space and the only degree of free- 
dom, which remains to be summed, are the path indices. 
We can get rid of these sums by defining vector currents 
j = (j + , — j_ ) and interpreting the 2-point cumulants as 
2x2 matrices C = Cp ± p 2 . This then allows to write the 
second-order term as j(C + JCcr 3 C)j* with the Pauli ma- 
trix a 3 taking into account the sign P. The subsequent Leg- 
endre transformation then leads to an iterative inversion of 
this coefficient yielding the inverse of the geometric series 
C ^„(«/<7 3 C)", which contains all hopping "chains". 
From the resulting Ginzburg-Landau functional 

T = tf Q [C~l - J^<5 Pq ,p,PJ % (5) 

we finally obtain equations of motion via extremization, 
since the currents in the original physical system vanish: 

Together with the complex conjugate of Eq. ©, there are 
in total four equations of motion, which simplify when the 
path-ordered quantities in T are rotated into the so-called 
Keldysh basis. Here one considers the linear combina- 
tions = (X + (t) + X_(t))/V2 and X A (t) = 
(X + (t) — X_(t))/\^2. Then the equation of motion, 
where the derivative is taken with respect to can be 
solved by assuming = 0, which implies Vl/ + = 
It turns out that the other equation of motion determining 
depends only on the retarded and advanced cumulants 
which are defined as the thermal average of multiple com- 
mutators times Heaviside step functions (cf. Ill5ll "). 
The non-trivial solutions of Eq. (O and its complex con- 
jugate can be divided into two classes which we work out 
at T = 0: 

1. Static solutions follow from the ansatz ^ S fc(w) = 
(5fc. ^(w) x I / oq . Eq. © then reduces to an algebraic equa- 
tion which determines \l/ eq . By inspecting the region in 
the parameter space, where \I/ eq starts to become non- 
vanishing, we find the phase boundary depicted for T = 
and 0<///{/<lon the left side of Fig. [TJ In the con- 
sidered first hopping order, this quantum phase diagram is 
identical to mean-field results [6] with a deviation from re- 
cent high-precision Monte Carlo data [19] of about 25%. 
Within a Landau expansion the second-order hopping con- 
tribution has recently been calculated analytically in Ref. 
fl 1 Oil decreasing the error down to less than 2%. A numeric 
evaluation of higher hopping orders has even been shown 
to converge to a quantum phase diagram which is indistin- 
guishable from the Monte-Carlo result 112.011 . 
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FIG. 1: Excitations in k — (1, 0, 0) direction are plotted in Al - CI for different values fj,/U and J/U, which are marked in the phase 
diagram (left). In the MI phase (green) and on the phase boundary (blue), the two T — modes can be interpreted as particle (dotted) 
and hole (dashed) excitations. At the tip of the Mott lobe (B), both modes become gapless, whereas for larger (smaller) chemical 
potentials /i only the gap of the particle (hole) mode vanishes. In the SF phase (red), the gapless mode turns into a sound mode, but 
also a gapped mode is present everywhere in the SF phase. The smooth transition from the MI excitation to the SF excitation is further 
analyzed in A 2 - C 2 and A 3 - C 3, where the effective mass m and the gap A of each mode are plotted as a function of J/U. The 
sound velocity c of the massless SF excitation, plotted in A3 - C3, vanishes at the phase boundary except at the tip, indicating the 
existence of a different universality class in this configuration. 



2. For obtaining dynamic solutions we perform the ansatz 
^sfc( w ) = <5fe,o^( w )^cq + S^i(uj) and consider only 
terms which are linear in S^f. With this we obtain alge- 
braic equations which determine the excitation spectra. As 
the corresponding analytic expressions are rather cumber- 
some, we restrict ourselves to discuss them graphically in 

Figiu 

In the MI phase, we find the gapped particle/hole excita- 
tions from mean-field theory [6]. At the phase boundary, 
we have to distinguish the tip from the rest of the lobe 
in accordance with the critical theory of the BH model 
fUEl. While at the tip of the mth lobe, i.e. at fj, m = 
■\frnljn + 1) — 1, both excitations become gapless and 
linear for small |fe|, off the tip at /i > ji m (/x < /i m ) 
only the particle (hole) mode becomes gapless and remains 
with a finite effective mass. In both cases the vanish- 
ing of the gap can be described by the exponential law 



A ~ [ J— Jp B (/u)] Ziy with the mean-field critical exponents 
v = 1/2 as well as z = 1/2 at the tip and z = 1 off the tip. 
Turning into the SF phase, the gapless mode rapidly loses 
its mass and has to be identified with the Goldstone mode 
which arises due to the broken U(l) symmetry. Indeed, 
within our Ginzburg-Landau theory it turns out for k — > 
and to — > 0, that the excitation 5^f stems from variations 
of the phase. Within a slave -boson approach it has even 
been shown in Ref. [7] that also for general wave vectors 
k phase variations dominate this excitation. This leads to 
density variations which make this mode sensitive to Bragg 
spectroscopy. Recently, the whole sound mode has been 
measured via Bragg spectroscopy far away from the phase 
boundary and could be well described via a Bogoliubov fit 
150. Surprisingly, also this regime deep in the SF phase 
turns out to be accessible with our Ginzburg-Landau the- 
ory. Expanding the Green's functions in U, in the lowest 



4 



non-trivial order they do not depend on temperature and 
lead to the equation of motion 



ih- 



dt 



J^Jij^j-^i-UVilVi] 2 , (7) 



which is the lattice version of the Gross-Piteevskii (GP) 
equation ll22ll . From this follows the Bogoliubov sound 
mode of a fully condensed system ]6J]. 
Additionally to that sound mode, however, also a gapped 
mode survives the quantum phase transition or arises again 
if we depart from the tricritical point. It can be smoothly 
mapped onto one of the respective MI modes, which is 
shown in the plot of both the effective masses and the gaps 
on the right side of Fig(T]The existence of such a SF gapped 
mode is in accordance with results obtained previously in 
Refs. flUI], however this finding has not yet been exper- 
imentally confirmed. The gapped mode seen via Bragg 
spectroscopy in Ref. [5|] can be well explained within a 
Bogoliubov-de Gennes ansatz considering second-band ex- 
citations due to finite temperature. In Ref. [7|] it is argued 
that the physical picture behind the gapped T = mode is 
an amplitude excitation which corresponds to an exchange 
between condensed and non-condensed particles at con- 
stant overall density. This reasoning is compatible with 
a numerical investigation of our full equations of motion 
(O where the gapped mode converges for small U to the 
constant dispersion u;(k) = 2/i. Furthermore, due to the 
absence of any density variation, this mode should be in- 
sensitive to Bragg spectroscopy. Although Eq. © does not 
allow pure amplitude excitations, we can back this inter- 
pretation of predominant amplitude excitations by observ- 
ing that the zero-momentum energy transfer at the phase 
boundary corresponds to the creation of a particle/hole pair. 
Via a modulation of the lattice potential, a finite energy ab- 
sorption at zero-momentum transfer has already been ob- 
served in Ref. J4J|. Although this might be seen as the first 
experimental signal for such a SF gapped mode, an unam- 
biguous quantitative identification is still lacking. 
Finally we compare our results with the ones obtained 
within a similar Ginzburg-Landau theory in imaginary time 
ifTTIl . Whereas at T = both formalisms yield identical re- 
sults, mismatches occur for finite temperature. This is sur- 
prising, since both formalisms are considered to be equiva- 
lent in equilibrium lfl~5l[l6ll . We conclude that the Keldysh 
formalism working with a purely real time-evolution con- 
tour is not able to produce the correct equilibrium config- 
uration of the full system for finite temperature. In a non- 
equilibrium system it might be justified to define the tem- 
perature only in the unperturbed initial state. For a system, 
however, that is supposed to relax into a new equilibrium 
state, a temperature change is expected, if the Hamiltonian 
describing the new equilibrium does not coincide the old 
one. From this reasoning follows that the imaginary part 
of the time-evolution contour must not be neglected II23I1 . 
Furthermore, the agreement of both formalisms at T = 



can be understood as a consequence of the Gell-Mann- 
Low theorem II24II stating that the systems remains in the 
ground-state, if a perturbation is adiabatically switched on. 
Since our ansatz has pushed this switching into the infinite 
past, no additional assumptions about its adiabatic proper- 
ties had to be made lfl6ll . 

The mismatch between our real-time calculation and the 
previous imaginary-time formalism suggests further inves- 
tigations where the time-evolution contour is extended to 
the imaginary time axis. Only then it is possible to deter- 
mine the temperature dependence of the excitation spectra 
which might serve as a thermometer for bosons in optical 
lattices UM- 
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